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Abstract 

The density of the hyperbohc metric on the complement of a rect- 
angular lattice is investigated. The question is related to conformal 
mapping of symmetric circular quadrilaterals with all zero angles. 
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0. Introduction. A famous theorem of Landau says that there exists L > 
with the property that for all functions / analytic in the unit disc, such that 
|/'(0)| = 1, the image contains a disc of radius L — e, for every e > 0. 

The largest constant L for which this is true is called the Landau constant, 
and its exact value is unknown. It is natural to conjecture that the extremal 
functions are universal coverings of complements to lattices in the plane. 

For lattices with the property that the largest disc in the complement has 
radius 1, one can consider the problem of maximizing |/'(0)| over all analytic 
functions mapping the unit disc to the complement of the lattice. Albert 
Baernstein II and J. Vinson [2] found that the hexagonal lattice and the 
universal covering with the property that /(O) is a center of the complement 
give a strict local maximum in this problem. 

This paper originated from an attempt of the author to prove a similar 
result for the class of rectangular lattices. A natural conjecture is that the 
square lattice must be extremal, and that for the extremal function /(O) is 
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a center of the complement. This is still unproved. However, the author 
feels that even the computation of the hyperbolic metric at the center of the 
complement of a rectangular lattice is an interesting problem. We will see 
that it is related to conformal mapping of certain symmetric circular quadri- 
laterals, whose all angles are zero, the problem which attracted attention of 
mathematicians and physicists during the whole XX century [5| IH 171 ISj ITU] . 

We consider a rectangular lattice {2muj + 2nuj'} where u = (lni^)/2, 
K > 1, and u' = ni. Let / be a universal cover of the complement of this 
lattice by the unit disc such that /(O) = u + u' is a center. We are interested 
in the quantity A{lnK) = |/'(0)| function of In i^. 

In our use of the standard notation of the theory of elliptic functions we 
follow [5] (see also [1]): r = u'/u, h = e^"; 9j{() is the j-th theta-function, 
9j = 9j{0), and Jacobi's Modular Function is k^. 

1. Denote 

fc=i^, a{k) = A{27rk) = Ai\nK) = \f'm. 

We may assume that / maps a circular quadrilateral Q (having zero angles, 
inscribed in the unit circle, symmetric with respect to reflections in coordi- 
nate axes) onto a fundamental rectangle R of the lattice, such that /(O) is 
the center of R. Then a simple symmetry and rescaling argument gives the 
functional equation 

a{k'^) = k-^a{k), k > 0. 

It is easy to see that in the limit when A; — )■ 0, / maps the unit disc onto the 
strip < \Qu\ < 27C, and we obtain 

a(0) =4. 

Together with the functional equation this implies 

a{k) ~ 4/c, /c — 7- oo. 
Differentiating the functional equation we obtain 

a'(l) = ia(l). 
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In the next two sections we find a(l) and a'(0). In section 4, we derive an 
"explicit expression" for a{k) for all k, and in section 5 discuss the compu- 
tational aspects. 

Finding the density of hyperbolic metric in the complement of a rect- 
angular lattice is equivalent to finding a conformal map from a rectangle 
onto a hyperbolic quadrilateral with zero angles. This is a classical problem 
investigated by Hilbert |1] and Klein [S], and in modern times in [7]. The 
mapping satisfies a Schwarz differential equation related to a Lame equation. 
The problem of finding this mapping for a given circular quadrilateral (not 
necessarily inscribed in a circle) requires determination of the so-called ac- 
cessory parameter which is a solution of a transcendental equation involving 
Hill's determinants. See [31 El ED] for results in this direction. The only paper 
I know where the accessory parameter was actually cornputed is [7] but it 
does not contain a rigorous justification of the algorithmic So we describe in 
sections 4 and 5 a convergent algorithm for the special case considered here. 



2. Finding a(l). We use the upper half-plane rather than the unit disc. 
Let T be the triangle in the upper half-plane with zero angles and vertices 
0, 1, oo. Quadrilateral Qi is the union of T with its reflection about its left 
vertical side. Now the "center" is at the point t = i, and we are looking for 
\f[{i)\, where 

fiir) = /(C), 

r(C) is a map prom the unit disc onto the upper halfplane, r(0) = i. Our 
function /i maps Qi onto a square Ri of side 27i. We split it into a compo- 
sition of two functions. 

fi = 9° 

where k"^ is the Modular Function of Jacobi. (In [T] this function has double 
notation, sometimes A, sometimes k^). maps T onto the upper half-plane 
sending (oo, 0, 1) to (0, 1, oo). Thus K,^{i) = 1/2. 

Our second component is a Schwarz-Christoffel map 

9{w) = c rc'^'ii-cr'^'dc 

Jo 

^The authors say on p. 217: "It should be emphasized that our remarks about the 
implicit equations are purely heuristic and that the actual computation proceeded, as it 
were, fortuitously without any a priori justification. 
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with 

27ry2 

C = , , , , ^ 1.981 , 
5(1/4, 1/4) 

where B is Euler's Beta-function. This g maps the upper half-plane onto 
the isosceles right triangle with legs 27i, and (7(1/2) is the middle of the 
hypotenuse. We have 

^'(1/2) = C-23/2 _ 3 gg7_ 
Using the notation from [H [5] , we have 

K' = ie2/es)\ 011,4, §5], 

03 = 1 + 2/1 + 2/1^ + 2/1^ + ..., (see [51 11,2, §6]), 



where 
Thus 



h = exp^nir). 



and dr/dC, = 2 at z = 0. 

For C, = T = h = e~'^ , Matlab gives 



|f(0)|=2|(«:7(z)|(7W2)^ 7.416. 



3. Finding a'(0). Let us make a preliminary map of the unit disc onto the 
horizontal strip \'^w\ < 2tt, so that our quadrilateral Q is mapped onto the 
quadrilateral Q2 bounded by two vertical segments of length 27r each, and two 
half-circles, orthogonal to the boundary of the strip, and zero corresponds to 
the center of Q2, which is also at ti^ = 0. We have 



dw 



dC 



4 at the point C = 0. 



Suppose that Q2 has width e > which is small. Then /2 maps Q2 into 
a rectangle R2, bounded by two vertical segments of length 2tt and two 
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horizontal segments, and also having small width. It is enough to estimate 
the width of this R2 because 

/^(O) « (width ofi?2)/e. 

Let us rescale both Q2 and R2 to the width 1, denoting the rescaled quadri- 
lateral by Q'2 and R2, and let G be the conformal map between them. Now 
we are interested in the ratio of the lengths of these quadrilaterals, which 
are large. If we cut Q2 and R2 by horizontal segments in the middle, the 
lower half of Q2 will be conformally equivalent to the lower half of R2 (as 
curvilinear quadrilaterals). Let us compare the restriction of G to these halts 
with the conformal map F of the triangle T (see section 1) onto a vertical 
half strip with vertices 0, 1, 00 (and right angles at and 1). Such map with 
the vertex correspondence (0, 1, 00) — >■ (0, 1, 00) is given by 

1 2 - k\t) 
F{t] — — arccos „ ' . 

Using the explicit expression for above and putting t — iy, we obtain 

F{iy)/i = y ln4 + o(l), y ^ +00. 

TT 

This means that i?2 is shorter than by 

-ln4 + o(l). 

TT 

(That G{iy) — F{iy) + o{l),y ^ 00 is quite evident, and this is easy to 
justify by using extremal length or some other argument). Passing back to 
the original R2 and Q2 we obtain that the width of R2 is 

e + — ln4 + o(e^). 

TT 

This means that a = |/'(0)| = 4|/^(0)| = 4 + AeTT-'^lnA + o(e), that is 

4 

a'(0) = — In 4 fa 0.5618. 



4. Representation of |/'(0)| = a{k) for arbitrary k. 
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Let p be the Weierstrass function of our rectangular lattice (with real 
period 2uj and pure imaginary period 2uj'). We set 

P{0 = \{p{C + uJ + uj')~e^). 

Then P is real on both real and imaginary axis, in fact it maps the rectangle 
-Ro = {0 < 3ftu < w, < < uj'/i} (one quarter of the fundamental 
rectangle) onto the lower half-plane (see Fig. 1) P is holomorphic in the 
closure of Rq, except at one point u + u' where it has a pole of second order. 




£3-62 
4 



£1—62 
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(0) 



C PiC) 

Figure 1. Conformal mapping by P. 



Consider the differential equation 



+ P{Ow = Xw, (1) 



where A = X{uj,u') is a real parameter to be specified later. 

Let A~ be the largest eigenvalue of ([T]) with boundary conditions w'{0) = 
w'{uj) = 0. Taking into account that P(C) < for C G [0,a;] (see Fig. 1), we 
conclude that 

A- < 0. (2) 

Putting ( = it in we obtain 

^ = (P(zt) - X)w. (3) 
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Let be the smallest eigenvalue of with the boundary conditions w'{0) = 
w'{u'/i) = 0. Using P{C) > for ^ G (0,0;'] (see Fig. 1), we obtain 

A+ > 0. (4) 

Let c and s be two solutions of ([I]) defined by 

/ c(0) s{0) 0] 
[ c'(0) s'iO) ) [O 1 )■ 

Then c{a),s{a) are real for a G [0,0;], while c{it) is real and s{it) is pure 
imaginary for t G [0,a;']. We introduce four real quantities: 

m = {s{uj)/c{uj) + s'{uj)/c'{u))/2, (5) 

r = is{co)/c{co)-s'iu)/c'{co))/2, (6) 

mi = -iisicu')/ciu') + s'icu')/c'{co'))/2, (7) 

n = -i{s{uj')/c{uj')-s'{uj')/c'{uj'))/2. (8) 
Proposition. There is a unique A G [A^,A^] such that the condition 
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ml + m? 

is satisfied. Then we have 



n (9) 



aik) ^ -^L= (10) 
Jm^ + nil 



Proof of the Proposition. Consider the Schwarz differential equation as- 
sociated with ([1]): 

2 

I =2(P-A). (11) 



F'" 3 F" 



It is well known (and easy to verify) that every solution of ( ITT]) is a ratio of 
two linearly independent solutions of ([T]). So we put F = s/c. 
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4.1. Then F is locally univalent, real on the real axis, and pure imaginary 
on the imaginary axis. 

4.2. We claim that for A G [A", A+], 

c(C)^0 when C e [0, U [0, w]. (12) 

This follows from the Sturm Comparison Theorem. 
Thus F is holomorphic on [0,0;] U [0,0;']. 

4.3. By the Symmetry principle, 

FiO = F{0 and F(-C) = -F(C). (13) 

It is well-known that F maps the right vertical side L' and the top horizontal 
side L of the rectangle Rq onto some arcs of circles on the Riemann sphere. 
We recall a simple proof of this. Consider, for example the right vertical 
side L'. As the right hand side of (fTTj) has period 2a;, and all solutions of 
(fTTj) are fractional linear transformations of each other, we conclude that 
F(C - 2u;) = A o F(C). Using L' - 2u = -77 and ((13]), we obtain 

-FXU) = Xo F{L'). 

Thus F{L') is fixed by an anticonformal involution, so it is an arc of a circle. 
Similar argument applies to F{L). 

4.4. It is known that the arcs F{L) and F{L') have exactly one common 
point zq and they are tangent at this point. We recall how this is proved. 
The coefficient P in ([1]) has one pole in the closure of Rq, namely at the point 
u + u', and the Laurent series at this pole has the form 

P(u; + u;' + C) = ^ + .... 

It follows that ([1]) has two linearly independent solutions near u + u' of the 
form 

^(w + w' + C) = C'/' + ... and V^(w + u;' + C) = C'/'lnC + .•• 
(see, for example, [6]). Thus F (which is a fractional-linear transformation 
of (p/ip) has a limit zq as ( ^ u + oj' from inside of Rq, and the angle between 
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F{L) and F{L') at zq is zero. As F{L) and F{L') are arcs of circles, Zq is the 
only point of intersection of these circles. 

4.5. We have shown in 4.1-4.4 that F maps the boundary dRo locally uni- 
valently onto a Jordan curve (consisting of a segment of the real axis, a 
segment of the imaginary axis, and two arcs of circles tangent to each other, 
one perpendicular to the real axis another to the imaginary axis). 








Figure 2. Conformal mapping by F. 



We claim that right circle has center at the point m and radius r, while 
the top circle has center irrii and radius ri (the definitions of these four 
numbers are given in ([5]). 

Let us verify the claim for the right circle. It is clear that the matrix 

s{u) s'{u) \ 
c{u) c'{u) J 
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is real. Introducing the pair of solutions (0, ip) of the equation 
by 

/ <P{u) {d<p/dt){u) \ / 1 \ 
I ^P{u) {d^/dt){uj) i I 1 i 



normalized 



where ( = it, 



we see that and ip are real on the vertical line {u + it : t & R} (as solutions 
of real differential equation with real initial conditions), and conclude 
that 



s 

c 



s{ijj)(t> + is'{(jj)ip, 
c{ijj)(j) + ic'{uj)ip. 



and thus the image F{L') belongs to the circle 



s[UJ) + s {u})ix 
ciuj) + c'{uj)ix 



: X e R 



whose center is at the point m and radius is r. 
The verification for the other circle is similar. 

Now we choose our parameter A in such a way that the common tangent 
line to the two circles passes through the origin. Elementary geometry shows 
that this happens iff the condition ([9]) is satisfied, see Figure 3. 



irrii 




Figure 3. Choosing the accessory parameter. 
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When A — A , the right circle becomes a vertical hne, and when A — t- A"*" 
the top circle becomes a horizontal line. It follows that the quantity e = 
m\/ \Jm? + nil ~ '"i changes sign when A changes on the interval [A~, A"*"]. It 
follows that there exists A such that condition dH]) is satisfied. In fact this 
value of A is unique but we do not use this additional information. 

For this value of A, F maps our rectangle i?o onto a quadrilateral Qq 
bounded by two straight segments F([0, a;]) C R, ^([0,0;']) C iR and two 
arcs of circle tangent at the point zq. Reflecting this quadrilateral Qq three 
times with respect to the coordinate axes we obtain a circular quadrilateral 
symmetric with respect to the coordinate axes, inscribed in the circle {z : 
\z\ = \zo\} and sides perpendicular to this circle. Thus F is the inverse to the 
universal covering of the complement of our lattice by the disc {z : \z\ < \zq\}, 
-F(O) = 0, F'{0) = 1, and it remains to verify that the radius \zq\ of this disc 
is given by the formula ffTOj) . This is clear from Figure 3. 

5. Remarks on computation. Our method of computation is based on 
the previous section. We are solving the Lame equation ([T]). We represent 
P in terms of theta-functions or in terms of elliptic functions of Jacobi (see 

my- 

^''^ 16u;2^i^i(C/(2u;)) 

where either product or series representations can be used for theta functions, 
they converge very well. Numerical solving of the equations ([T]) and on 
the real line involves only computations with real numbers. 

As Matlab does not have standard routines for theta functions, we can 
use an expression of p in terms of Jacobi elliptic functions. Matlab uses 
AGM to compute them which is probably as effective as theta functions. 

To find the accessory parameter A we use a "shooting method", dissecting 
the interval [A~,A+] dyadically. 

The result is that a(l) ~ 1.4163 and the graph of |/'(0)| for the rect- 
angular lattice with sides 2a;, 2a;', 4a;^ + 4a;'^ = 1, as a function of 2a; is 
this: 
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1.7 



1.65 



1.6 



1.55 



1.5 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



We see a maximum at the point 1/a/2 and this maximum equals 1.6693. 
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